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Abstract 

Using a method previously developed, based on the Mellin-Barnes transform, 
we reconstruct the two-point correlators in the vector, axial, scalar and pseu- 
doscalar channels from the Taylor expansion at q 2 = 0, the threshold expansion 
at q 2 = Am 2 and the OPE at q 2 — > — oo, where m is the heavy quark mass. The 
reconstruction is analytic and systematic and is controlled by an error function 
which becomes smaller as more terms in those expansions are known. 



1 Introduction 



Vector, axial-vector, scalar and pseudoscalar two-point correlators of heavy quarks are 
very instrumental functions for the determination of fundamental QCD parameters 
such as the strong coupling constant, a s , and the heavy quark masses, m. Presently, the 
input for these determinations comes both from the fact that QCD has been "solved" 
in e + e~ experiments [U 121 El H] and also from the fact that QCD may be "solved" nu- 
merically on the lattice [5]. Comparison between these two types of approaches results, 
of course, in a nontrivial test of the theoretical ideas involved and serves to assess the 
progress in our detailed understanding of QCD. 

The fact that the quark is heavy brings about the welcome simplification that 
perturbation theory is a valid approximation. However, even in this case, and already 
at 0(a 2 ), present state-of-the-art perturbative calculations do not allow a complete 
knowledge of these two-point functions for all values of the momentum q 2 . Instead of 
this, only expansions around specific values of q 2 have been possible to obtain. In Refs. 
[6]- [20] we summarize the vast literature showing the titanic effort of several groups 
to compute the expansion of these two-point functions at low energy (q 2 = 0), at the 
production threshold (q 2 = 4m 2 ) and at high energy (— q 2 — > oo). In some particular 
example, even up to 30 terms have been calculated. Such is the case of the low-energy 
expansion of the vector two-point function at 0(a 2 ) in Ref. [El [9]. However, more 
often than not, and certainly at 0(a^), we are limited by the fact that just a few terms 
in each expansion have actually been computed. For 0(a^), a nice summary of all the 
information about these three expansions, gathered up to date, can be found in Refs. 

[2D 122!. 

Given the limited amount of information, a full reconstruction of the corresponding 
two-point function can only be approximate. In the pioneering work in Refs. [23]-|27j. 
Pade Approximants [281 29J were proposed as a method to simultaneously resum the 
three individual expansions. However, a Pade Approximant is a rational function of q 2 
and, therefore, cannot reproduce either the logarithms of the OPE or the square roots 
of the threshold expansion. Consequently, one must first subtract from the original 
Green's function all this nonanalytic behavior with the help of some simple, but physi- 
cally motivated, functions. Even so, some of the Pades constructed turn out to contain 
unphysical features, and should be discarded by means of some type of prescription. 
This results in some degree of arbitrariness inherent to the method and, of course, in 
a corresponding systematic error, as was pointed out in, e.g. Ref. [21J. A common 
practice then is to estimate this systematic error by some variation over the different 
more-or-less arbitrary choices made along the process of construction of the Pades. 

Of course, whenever possible, checks have been made to assess the reliability of this 
Pade-inspired method. Although normally there seems to be a (perhaps surprisingly) 
high level of success, we think it is important that the results may be compared with 
those obtained by a completely different technique, such as the one which will be used 
in this work. As a matter of fact, Pades are known to be notoriously unpredictable ap- 
proximants, with seemingly capricious properties of convergence. A notable exception 
to this rule is when the function to be approximated is a Stieltjes function, i.e one with 
a positive definite spectral function, for which there is a theorem assuring convergence 
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of certain Pade sequences [28], [29] . Notice however that, although the complete vacuum 
polarization functions (12. 2p are Stieltjes, a fixed order in the a s expansion need not be. 
In particular, as one can see in Figs. [3] below, the 0(oP s ) spectral functions are clearly 
not positive definite. This means that the convergence of the Pade-inspired method is 
not guaranteed. As a matter of fact, past examples in similar contexts do exist where 
the method produced certain results with an expected accuracy, and a comparison to 
the exact answer revealed that this accuracy had been overestimated [50] . 

In Ref. [21] an alternative method has been presented. Unlike the Pade method, 
which is heavily numerical, the method in [31] is completely analytic. The approxima- 
tion is also systematic (e.g. it obeys a simple counting rule) and expresses the expansion 
of any vacuum polarization function as a unique combination of polylogarithms plus 
a known polynomial in the conformal variable u. This variable maps the cut complex 
plane in q 2 onto a unit disc. Furthermore, a parametrization of our ignorance about 
the first term not included in the approximation allows one to make an estimate of the 
systematic error incurred. 

In this work, we will obtain results for the vector, axial-vector, scalar and pseu- 
doscalar vacuum polarization functions at 0(af) and make a detailed comparison 
with the corresponding results using the Pade-inspired method, when available. 



2 Definitions 



We define the generalized currents in the four channels: scalar (s), pseudoscalar (p), 
vector (v) and axial-vector (a) as 



and the corresponding vacuum polarization functions as 



(2.1) 



q 2 IP'V 



q 2 g, v ) W^{q 2 ) + q,q u YiT{q 2 ) 



(fx e igx (0|T 
d 4 x e iqx (0|T 



~s(x)s(0) 
p(x)p(0) 

a M (a:)a I/ (0) 



10) 



|0>. 



(2.2) 

We will refer to this set of four polarization functions globally as IP where the index 
X will run over the four channels, i.e. x = s iP) v i a - 

The longitudinal polarization function U v £ a vanishes in the vector channel and, in 
the axial channel, is related to the pseudoscalar function through a Ward identity, so 
it will not be considered any further. Moreover, topologies in which massless particles 
occur as intermediate states have to be treated separately since they give rise to dif- 
ferent analytic properties of the polarization function in the complex plane. Leaving 
these contributions with massless cuts for future work, here we will concentrate on 
those perturbative Feynman diagrams which cause the polarization functions IP to 
have a cut in the complex plane starting at q 2 = 4m 2 , where m is the (pole) mass of 
the heavy quark considered. 



1 Some of these results for the vector channel were already obtained in Ref. [31] 
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Figure 1: Conformal mapping between z and cu, Eq. (I2.3j) . 



Let us define the new variables z = q 2 /Am 2 and w as 

4w 1 - v/T^ 



(2.3) 



This change of variables maps the cut z plane into a unit disc in the u plane, as we can 
see on Figure [U The physical cut z G [1, oo[ is transformed into the circle \ou\ — 1 . The 
points z = into lo — 0, z — 1 into u — 1, the limit z — >■ +oo ± into a; — > — 1 ± ie, 
and z — oo into u — >■ — 1. 

In terms of the variable z, the polarization function II x (z) has the following prop- 
erties: 

1. It admits a convergent expansion when z — > since it is analytic in the disc 
Izl < 1: 



U x (z) 



\z\<l ^ 



n)z 



(2.4) 



71=0 



2. It has a cut along 1 < Re(z) < oo. In the limit 2 — >■ 1, it can be expanded 
(threshold expansion) in the form 



z— s>l ^ — ' 

p,k 



(2.5) 



3. It admits and expansion (Operator Product Expansion) for z — > —00 such as 

IF(z) ~ ^^(p,A:) llog fc (-4z), (2.6) 



p.k 



On the other hand, in terms of the variable u the polarization function II x (u;) = 
Il x (4a;/(1 + u) 2 ) obeys the following three expansions: 
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1. The Taylor expansion (u — > 0): 

oo 

IF( W ) = V(F(n) w n . (2.7) 
M<i ^ 

n=0 

This expansion has radius of convergence \u>\ < 1. 

2. The threshold expansion (u — > +1) and the Operator Product Expansion (OPE) 
(w — > — 1) can be combined into: 

8*( W ) ~ V^W(A,p)(lTo;) A log p (lTa;) , (2.8) 

A,p 

where the index p is a (positive) integer but the index A may be integer (OPE), 
and also half-integer (threshold). 

The Taylor expansions of the function U x (w), Eq. (12. 7p . and that of the function 
n x (z), Eq. ( 12.41) . are related. The relation between the two is given by 

t^r(2p)r(n + i- P ) w ' 1 ; 



r (n) r(| + n) 



c*(n) = V ' rF( P ) . (2.10) 

v ; ^r(i + n- P )r(i + n + p) w v ; 

It was shown in Ref. [21] that the asymptotic behavior of the Q x (n) for large n 
is solely determined by the coefficients fl x ^(X,p) in Eq. (12.81) encoding non-analytic 
behavior, i.e. when A < 0,Vp, or when A > 0,p ^ 0. Furthermore, this asymptotic 
behavior is very precocious, yielding accurate results already for n > 2, 3. The com- 
bination of these two features allows the reconstruction of the vacuum polarization 
function by carrying out the sum in Eq. (12.71) using the exact value for the first few 
coefficients Vt x {n) up to a given value, say n = N*, and then using the asymptotic 
expansion for n > N*. 

More specifically, the coefficient Q x admits the following expansion 



p,fe 



where the coefficients and (3* k are in one-to-one correspondence with the coef- 
ficients Q x( - + \\,p) and Q x (~\\,p) (respectively) in a systematic way. The precise 
connection is best obtained through: i) the Mellin-Barnes representation 



c+ioo 



n*(± e -*) = o / — r s r( s ) £ (±!) n fiX («) > ( 2 - 12 ) 

c-ioo n=1 

where c is a constant located in the region of analyticity of the integrand in the complex 
s plane, also known as the "fundamental strip" ; and ii) the Converse Mapping theorem 
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dictionary, by which a function, f(t), and its corresponding Mellin transform, Ai[f](s), 
are related through 

c+ioo 

f°° f ds 

M[f}(s) = j dtt^fit) ^ f(t)= J —r s M[f](s) (2.13) 

c—ioo 

where the symbol x means that one sums over all the negative powers of the Laurent 
expansions of the function around every pole (i.e. the "principal part" of the mero- 
morphic function). Feeding Eq. f |2. 1 2[) with the expansion (12. lip and matching onto 
the expansion for n x (±e -i ) on the lefthand side as given by (12. 8p with the help of 
the dictionary (I2.13|2.14p . one obtains all the coefficients a v ^ and f3 Py k in terms of the 
threshold and OPE coefficients A x (p, k) and B x (p, k) in Eqs. (I2.5|2.6p in a completely 
systematic way, order by order in the expansions. For instance, if the leading term in 
the threshold expansion is given by the coefficient A(— 1,0) in Eq. (I2.5P (as indeed it 
happens for \ — v,p, see the third row in Tabled]), one then obtains that the leading 
asymptotic behavior is given by ( 12. lip with a-1,0 = 4 A(— 1,0) (see the correspond- 
ing entry in Table [3]). The relation between the coefficients A and B in the OPE 
and threshold expansion (12.5112.61) and the coefficients a, in the asymptotic expansion 
(12. lip is exact. Although we will give numerical results for the sake of brevity, it is 
clear that the whole approach is fully analytic. 

3 Results 

All the relevant expansions have been conveniently listed in Refs. [211 122] for the 
four polarization functions Il x with x = v,a,p,s. Their corresponding coefficients 
for the threshold expansion (12. 5p and the OPE (12.61) can be extracted from these 
references and we collect them in Table [T] (where rii stands for the number of light 
flavors). Alternatively, we can also express the four polarization functions in terms of 
the variable u in Eq. (12. 3p . After expanding, the corresponding coefficients for the 
OPE and threshold expansion (12.81) can also be extracted and we collect them in Table 

m 

Using the matching condition (12.121) . one can deduce the values for the coefficients 
a* k and f3* k in the asymptotic expansion (12.111) . We list them in Table [3J This 
asymptotic expansion can then be compared to the exact values of the first coefficients 
Q x (n) in Eq. (12. 9p . where the C x (n) are known exactly by means of diagrammatic 
methods. This is done in Figures [21 As one can see, the asymptotic expansion starts 
being a good approximation already for n > 2, 3. Note, in particular, the pseudoscalar 
channel where this comparison may be done up to n = 4. 

Once the coefficients a Pi k and f3 P} k are known, one may reconstruct the vacuum 
polarization function Il x (u;) by directly carrying out the sum in Eq. (12.71) . splitting it 
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IF 


fiz = 3 


m = 4 


ni = 3 


ni = 4 


ni = 3 


ni = 4 


ni — 3 


ni = 4 


A(-l 0) 


2.6364 


2.6364 


o 


o 


2.6364 


2.6364 










-25 2331 


-24 5549 


o 


o 


-22 9364 


-22 2582 










-7.7516 


-7.1774 


o 


o 


-7.7516 


-7.1774 










-11.0654 


-9.1783 


-0.7311 


-0.7311 


-3 8928 


-1.8142 


-1.0966 


-1.0966 




1.4283 


1 4732 


n 

U 


n 

u 


3 9655 


3 2826 


n 

U 


n 

U 


A(0, 3) 


-0.4219 


-0 361 7 


n 

U 


n 

u 


-0.4219 


-0 361 7 

U. 0\J i- t 


n 

U 


n 

U 


5(0, 1) 


n nfiQQ 

u.uuyy 


U.U± I o 


OfiQQ 


01 7^ 

U.U± 1 O 




1 ^1 94 




1 ^1 94 


5(0, 2) 


1911 


OQQO 


1911 


OQQO 




O. ( i 1^ 


4 ^OQ 


q 77I O 


-8(0, 3) 


u.uouu 


u.UO 1VJ 


O^fifi 
u.uouu 


031 


1 1 841 




1 1R41 


1 Ofifi& 

-X .UUUO 


B(0 4) 


o 


o 


o 


o 


0.0797 


0.0723 


0.0797 


0.0723 


5(1,1) 


-3.7567 


-3.1062 


-3.3602 


-2.7103 


-12.0583 


-10.3098 


-18.0487 


-15.3645 


5(1,2) 


2.1173 


1.8983 


-2.5325 


-2.1623 


3.0251 


2.8886 


-5.5011 


-4.6759 


5(1,3) 


-0.3189 


-0.2894 


0.8652 


0.7774 


0.4023 


0.3532 


3.7238 


3.4077 


5(1,4) 








-0.0797 


-0.0723 


-0.1573 


-0.1467 


-0.4719 


-0.4402 


5(2,1) 


-5.1301 


-4.3252 


4.5072 


3.7977 


-5.8848 


-5.3034 


11.8587 


10.3358 


5(2,2) 


0.3182 


0.3487 


-0.5490 


-0.5518 


-1.1566 


-0.9612 


-4.5366 


-4.3615 


5(2,3) 


0.4015 


0.3627 


-0.3190 


-0.2806 


1.3083 


1.2114 


-1.1407 


-1.0030 


5(2,4) 


-0.0787 


-0.0734 


0.0787 


0.0734 


-0.1565 


-0.1475 


0.4695 


0.4426 



Table 1: Threshold expansion and OPE coefficients in Eqs. (12. 5112. 6ft for ni = 3,4; 
where rii stands for the number of light flavors. 
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ni = 3 


ni = 4 


ni = 3 


ni = 4 


m = 3 


ni = 4 


ni = 3 


n/ = 4 




-2,0) 


10.5456 


10.5456 








10.5456 


10.5456 








fi*(+) 


-1,0) 


18.4286 


18.6642 








13.8352 


14.0708 








fi*(+) 


-bl) 


31.0063 


28.7095 








31.0063 


28.7095 










(0,1) 


-50.4189 


-45.0511 


-1.4622 


-1.4622 


-46.2611 


-40.3563 


2.1932 


2.1932 


n*(+) 


(0,2) 


12.7314 


11.9097 








20.0801 


19.1472 








n*(+) 


(0,3) 


-3.3750 


-2.8935 








-3.3750 


-2.8935 








nx(-) 


(0,1) 


0.4872 


0.2979 


0.4872 


0.2979 


-10.2951 


-7.9785 


-10.2951 


-7.9785 


n*c-> 


(0,2) 


-0.7349 


-0.6359 


-0.7349 


-0.6359 


-7.3644 


-7.0625 


-7.3644 


-7.0625 




(0,3) 


0.2932 


0.2481 


0.2932 


0.2481 


2.3993 


2.1162 


2.3993 


2.1162 




(0,4) 














1.2757 


1.1575 


1.2757 


1.1575 




(2,1) 


0.1671 


0.2416 


-2.2712 


-1.6012 


-1.5908 


-1.2726 


-3.3365 


-2.0631 


fF(-) 


(2,2) 


0.7552 


0.6947 


-0.7669 


-0.7815 


4.5968 


4.2651 


0.0063 


-0.0398 


fi*(~) 


(2,3) 


-0.6378 


-0.5787 


-0.0380 


-0.0497 


-1.4088 


-1.3911 


-1.7440 


-1.7921 


fix(-) 


(2,4) 








0.3189 


0.2894 


0.6292 


0.5870 


1.8876 


1.7610 




(4,1) 


0.7588 


0.7491 


-3.7157 


-3.0118 


-0.1969 


0.1395 


-8.5663 


-6.9123 


nx(-) 


(4,2) 


-0.7823 


-0.7052 


-0.3778 


-0.4278 


-0.9318 


-0.9492 


-3.1081 


-3.0457 


n*c-) 


(4,3) 


-0.4024 


-0.3532 


0.6421 


0.5518 


0.0547 


0.0036 


1.2480 


0.9595 


fi x <-) 


(4,4) 


-0.0787 


-0.0734 


0.3976 


0.3627 


0.4727 


0.4395 


2.3571 


2.2036 



Table 2: Threshold expansion and OPE coefficients in Eqs. (12. 8p for ni = 3,4 (number 
of light flavors) . 
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Figure 2: Comparison of the asymptotic expression VL\ s {n) (blue crosses) and the exact 
values Vt^ s {n) (red circles), for x — s ,P, v , a , extracted from Ref. [2lT [22] . 
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ni = 4 


ni = 3 


ni = 4 


m = 3 


m = 4 


" X i,o 


10.5456 


10.5456 








10.5456 


10.5456 








"0,0 


-11.0769 


-12.6382 








-6.4835 


-8.0448 








"0,1 


31.0063 


28.7095 








31.0063 


28.7095 








a l,0 


36.3318 


33.0585 


1.4622 


1.4622 


40.6575 


36.7189 


2.1932 


2.1932 




37.1514 


33.8404 








51.8488 


48.3155 








1,2 


10.1250 


8.6805 








10.1250 


8.6805 










-0.1819 


-0.0555 


-0.1819 


-0.0555 


9.9493 


6.9861 


9.9493 


6.9861 




-2.4852 


-2.1312 


-2.4852 


-2.1312 


-43.1187 


-39.6735 


-43.1186 


-39.6735 




-0.8795 


-0.7444 


-0.8795 


-0.7444 


1.6381 


1.6688 


1.6381 


1.6688 
















5.1027 


4.6298 


5.1027 


4.6298 


/3 3 X ,o 


-10.4385 


-9.7282 


26.2458 


22.9826 


3.1298 


1.1687 


93.7790 


83.0590 


Ph 


-4.7750 


-4.2501 


-19.8617 


-18.4878 


-53.4944 


-50.0465 


-137.2835 


-129.2810 




3.8270 


3.4724 


-6.8349 


-6.1103 


0.8960 


1.1335 


-24.9630 


-22.4605 


/3 3 x 3 








2.5513 


2.3149 


5.0337 


4.6960 


15.1011 


14.0879 


/f,o 


-70.9277 


-63.8573 


100.2171 


89.1103 


-115.8498 


-108.3750 


440.1394 


399.7520 


/f,i 


56.3093 


53.6862 


-72.4918 


-68.9185 


62.1988 


60.2675 


-512.9781 


-487.4843 


Ph 


20.9951 


19.0619 


-29.3263 


-26.2676 


38.4395 


35.8466 


-129.9058 


-118.3556 


Ph 


-7.5506 


-7.0439 


10.1019 


9.3589 


-9.9903 


-9.4668 


60.1732 


56.5763 



Table 3: Values for the coefficients in the asymptotic expansion (12.111) for rii = 3,4. 
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Table 4: Error functions. 



in three different terms: 

oo N* 

n*(w) = £ fl^n) + X] - w n + £ X {N\ u) , (3.1) 

n=l n=l 

where 

oo 

e*(N*,u,)= Yl [W(n)-nl s (n)]u n . (3.2) 

n=N*+l 

Notice that the first N* coefficients, which are known from the Taylor expansion (12. 7p . 
have been included exactly, using the asymptotic expansion Q^ s (n) only from the value 
N* + 1 onwards. The function S x is an "error function" that will have to be estimated. 
As the expression (13.21) clearly shows, S x (N*,u) becomes smaller as N* grows, for any 
| Co? | < f. One may also expect the error function to decrease as more terms in the 
series (12. If p are included. This means that our approximation will become better as 
more terms from the three expansions (I2.4ti2.6l) are known, as we explicitly show in the 
Appendix for the pseudoscalar channel. 

In this way, we could obtain the following approximation to the vector vacuum 
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polarization function [3T] : 

ff(u;) = a_i 7z - X9 + «o,o t~~ a ,iLi (1) (0,w) - a 1(Q log(l - co) 

(1 — tu) z 1 — UJ 

- a 1A Li«(l, u) + a 1>2 U^(l, u) - P lj0 log(l + u>) - /? M Li«(l, -u) 
+ ft , 2 Li< 2 > (1, -u) + P 3 , Li(3, -w) + /3 3 ,2 Li (2) (3, -w) + /3 5 ,o Li(5, -w) 

- /? 5 ,i Li« (5, -w) + /3 5 ,2 Li (2) (5, -w) - A;,3 Li^ (5, -w) 

AT* 

+ J] -fi A5 (n)] w n + ^(A^*,w) (3.3) 

n=l 

in terms of the polylogarithmic function 

IA( 8 ,u) = -P*— f'-log^f-) tog(l-rra;), (3.4) 



r(s — i) Jo x '" \s 

and its derivatives 



Li^( S ,c) ee ^Li( S ,c) = (-l) P f:^ • (3-5) 



n=l 



Analogous expressions can also be easily obtained for the other 3 channels: it is a simple 
matter to associate the coefficients a Pt k and f3 Pj k in Table [3] with the corresponding 
polylogarithm, after the pattern seen in Eq. (13. 3p . 

For a given number of known Taylor coefficients, N*, the error function (13. 2p may 
be obtained after a conservative estimate is made of the expected behavior for the 
difference [fl x — Q\ s ] in Eq. (13. 2D . This difference is expected to be of the order of the 
first unknown term in the asymptotic expansion (12.111) . Since our present ignorance 
is dominated by the a x coefficients in this expansion, we choose to parameterize this 
ignorance in the form of an effective a x . This effective coefficient is shown in Table 
H] for all the channels. Since these effective coefficients affect the value of the non- 
logarithmic terms of the OPE, i.e. terms of the form B x (p, 0) in Eq. ( 12. 6p . we have 
checked that these estimates are not in conflict with the value of those. 

In order to reduce the error even further, one would need to increase the value of 
N*, i.e. the number of Taylor coefficients in (12.41) . or the number of nonanalytic OPE 
coefficients B x (p, k ^ 0) in Eq. (12. 6p but, above all, the most efficient way would be to 
compute more nonanalytic terms of the threshold expansion, i.e. A x (p,k ^ 0) in Eq. 
(12. 5p (see the Appendix for some illustration on this point). 

Having an explicit reconstruction of the vacuum polarization function, e.g. (13. 3p . 
and an estimate of the error function, e.g. Table HI it is now a simple matter to also 
reconstruct the spectral functions through, for instance, the following representation: 

im u(s ' e ' v)= w) S? x ^iC^U • - if - % ■ (3 - 6) 

We show in Figs. |3] the result for the spectral functions in the four channels, as a 
function of the quark velocity v = \J\ — 1/z. The error on each of these curves 
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Figure 3: Imaginary part of the vector, axial-vector, scalar and pseudo-scalar polar- 
izations as a functions of the velocity v = a/1 — TJz. We also show the threshold 
expansion (in dashed green) at low v, and the OPE (in dashed black) at v near unity. 



is roughly of the order of the width of the curve itself. These results are in rather 
good qualitative agreement with the spectral functions found with the Pade method 
in Refs. [2T] and [22], as can be judged from the corresponding plots. A more detailed 
comparison would require knowledge of the error function associated with the Pade 
method, which has become available to us only in the vector channel^] For instance, 
at v = 0.75, we obtain 12vrulmn^ ) (l/l - v 2 ) = -246.7(8.8), whereas Ref. [21] obtains 
—219.0(3.4), and Ref. [22] gets —214.6, with an error that we have not been able to 
quantify in any detail. 

Knowledge of the asymptotic expression (12. lip allows us to predict the Taylor 
coefficients C x (n) through Eq. (I2.10p by separating the first iV* known terms in this 
sum and using the asymptotic expression from the term iV* + 1 onwards. Since these 
coefficients are of much phenomenological interest for the present determination of the 
heavy quark masses pQ-[5] we also list them in Tables [MS together with the analogous 
results found in Ref. [22] (see also Ref. [21] )• Our results confirm those of [22], albeit 
the errors quoted in [22] are always significantly smaller. 

2 A. H. Hoang and V. Mateu, private communication. 
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1 a 2 

C™(n) 


This work 




Ref. 


22 




n l — 


3 


Til — 


4 


n\ = 


3 


Til 


4 


n 




Error 




Error 




Error 




Error 


1 


366.1748 





308.0188 





366.1 / 48 





moo 
308.U188 





o 
z 


ool . OUy 1 


u 


ooU.OooO 


u 


ool.oUyl 


U 


OOU.OOOO 


U 


3 


385.2331 





338.7065 





385 2331 


o 


338 7065 


o 


4 


382.7 


0.5 


339.7 


0.5 


383.073 


0.011 


339.913 


0.010 


5 


378.0 


1.2 


337.7 


1.2 


378.688 


0.032 


338.233 


0.032 


6 


372.5 


1.8 


334.5 


1.8 


373.536 


0.061 


335.320 


0.063 


7 


367.0 


2.3 


330.9 


2.3 


368.23 


0.09 


331.90 


0.10 


8 


361.5 


2.7 


327.2 


2.7 


363.03 


0.13 


328.33 


0.14 


9 


356.4 


3.1 


323.5 


3.1 


358.06 


0.17 


324.78 


0.18 


10 


351.6 


3.4 


320.0 


3.4 


353.35 


0.20 


321.31 


0.22 



Table 5: Vector channel. 



^ C a (n) 


This work 




Ref. 


22 




ni = 


3 


«! = 


4 


«; = 


3 


n-l = 


4 


n 




Error 




Error 




Error 




Error 


1 


165.1328 





138.1938 





165.1328 





138.1938 





2 


105.1185 





90.0956 





105.1185 





90.0956 





3 


75.5564 





65.5198 





75.5564 





65.5198 





4 


57.7 


0.03 


50.42 


0.03 


57.7298 


0.0029 


50.4287 


0.0042 


5 


46 


0.1 


40.4 


0.1 


46.005 


0.009 


40.397 


0.013 


6 


37.8 


0.1 


33.3 


0.1 


37.813 


0.017 


33.338 


0.024 


7 


31.8 


0.2 


28.1 


0.2 


31.825 


0.025 


28.151 


0.036 


8 


27.3 


0.2 


24.2 


0.2 


27.291 


0.034 


24.206 


0.048 


9 


23.7 


0.2 


21.1 


0.2 


23.759 


0.042 


21.123 


0.059 


10 


20.9 


0.2 


18.6 


0.2 


20.943 


0.049 


18.658 


0.069 



Table 6: Axial channel. 



^ CP(n) 


This work 




Ref. 
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«i = 


3 


ni = 


4 


ni = 


3 


ni = 


4 


n 




Error 




Error 




Error 




Error 


1 


16.0615 





8.6753 





16.0615 





8.6753 





2 


230.9502 





199.8289 





230.9502 





199.8289 





3 


320.5093 





283.8922 





320.5093 





283.8922 





4 


359.1116 





321.5253 





359.1116 





321.5253 





5 


376.3 


0.3 


338.2 


0.3 


376.3673 


0.0023 


339.2386 


0.0021 


6 


383.4 


0.9 


345.9 


0.9 


383.6206 


0.0084 


347.4338 


0.0075 


7 


385.4 


1.5 


348.8 


1.7 


385.794 


0.018 


350.695 


0.017 


8 


385.7 


2.2 


349.1 


2.4 


385.250 


0.032 


351.252 


0.029 


9 


382.5 


2.8 


347.9 


3.1 


383.215 


0.048 


350.278 


0.044 


10 


379.5 


3.4 


345.9 


3.8 


380.360 


0.066 


348.424 


0.061 



Table 7: Pseudoscalar channel 



^ C s (n) 
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[22J 




ni = 
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ni = 


4 


ni = 


3 


ni = 


4 


n 




Error 




Error 




Error 




Error 


1 


-2.0665 





-3.9663 





-2.0665 





-3.9663 





2 


59.9301 





49.7941 





59.9301 





49.7941 





3 


69.5687 





59.9811 





69.5687 





59.9811 





4 


64.8 


0.3 


56.6 


0.3 


64.641 


0.014 


56.534 


0.014 


5 


57.4 


0.8 


50.6 


0.8 


57.168 


0.043 


50.399 


0.041 


6 


50.4 


1.2 


44.7 


1.2 


50.069 


0.081 


44.374 


0.076 


7 


44.4 


1.5 


39.5 


1.5 


43.95 


0.12 


39.10 


0.12 


8 


39.3 


1.8 


35.1 


1.8 


38.81 


0.16 


34.64 


0.15 


9 


35.1 


2.1 


31.3 


2.1 


34.52 


0.20 


30.89 


0.19 


10 


31.5 


2.3 


28.2 


2.3 


30.93 


0.24 


27.73 


0.22 



Table 8: Scalar channel 
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4 Comparison with the Pade method 



By construction, the asymptotic expansion presented in this work in Eq. ( 12. lip contains 
a finite sum which stops at a given order in k and p. These indices are determined 
by the number of nonanalytic coefficients, AOp, k) and B(p, k), which are known from 
the OPE and threshold expansions ( I2.5ll2.6p □. This means, in particular, that there 
is no way to predict the next term in the OPE or the threshold expansion from the 
previously calculated terms. This ignorance must be completely encoded in the error 
function 13.21 (see Table H]). The situation is similar to what happens in a generic Taylor 
expansion, where knowledge of n terms does not give any information about the term 
n + 1. 

This is in contradistinction to the Pade method. In this method, after all the nec- 
essary choices are made, one obtains a family of functions which, upon reexpansion 
around z — > — oo and z — > 1, yields a prediction for an infinite sequence of coefficients 
in both the OPE and the threshold expansion. Therefore, this means that the func- 
tions reconstructed via the Pade-inspired method contain more information than those 
reconstructed via the method of the present work. This could be consistent with the 
fact that the errors quoted within the Pade method always come out to be smaller, as 
shown in the Tables [5HH1 An important question to answer, however, is whether the 
prediction for the OPE and threshold expansion coefficients made by the Pade method 
may be considered sufficiently accurate. To clearly answer this question, the best will 
be to compare to the exact result, once the exact result for the next coefficient becomes 
available using standard diagrammatic techniques. 

In order to facilitate this comparison in the future, we use the one example of 
reconstructed function provided in Ref. |22j for each channel, and collect the prediction 
for the first coefficients in both the OPE and threshold expansion. This we do in 
both the z and u variables in Table |9] and Table HOI respectively. We also deduce, 
through the matching condition ( 12.121) . the corresponding a p ^ and /3 Pt k coefficients in 
the asymptotic expansion (12. lip , and list them in Table [TT1 Finally, we repeat the same 
exercise but using the family of reconstructed functions provided by the authors of Ref. 
[2T] . In this reference they considered the vector channel only, but they produced a set 
of Pade-reconstructed functions which passed all the good-quality criteria imposed by 
these authors. The results for this case are listed in Table O The two numbers quoted 
in every entry of this table correspond to the maximum and minimum value obtained 
within the set of reconstructed functions. When the two numbers are identical, this 
means that the number predicted turns out to be unique. By comparing the results 
shown in these tables one can see that sometimes there is complete agreement between 
the two analyses (like in the case of the A(l/2,2) coefficient), but other times the 
agreement looks more questionable (like in the case of the B(3, 4) coefficient). It will 
be very interesting to compare these predictions to the exact results, once they will 
become available in the future. 

3 By nonanalytic coefficients, we mean those A(p, k) and B(p, k), with p < 0, Vfc or p > 0, k ^ 0, 
i.e. those which accompany non-analytic behavior in 1 — z and 1/z, respectively. 
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a 


V 


s 




n ; = 3 


m = 4 


71; = 3 


ni = 4 


n\ = 3 


m = 4 


m = 3 


n t = 4 


A(l/2,0) 


-143.0996 


-76.2099 


-48.2359 


-44.0331 


-156.1437 


-161.1861 


-78.7753 


-70.4228 


A(l/2,1) 


-13.7341 


-6.4188 








-15.8236 


-15.6340 








A(l/2,2) 


3.2229 


2.7631 








3.2229 


2.7631 








B(3,l) 


-1.2191 


-1.2247 


0.9647 


0.8049 


0.4098 


0.2578 


4.3843 


3.8451 


B(3,2) 


-0.2746 


-0.4585 


0.1063 


0.0872 


-1.5633 


-1.4233 


0.5702 


0.4616 


B(3,3) 


0.2560 


0.4425 


-0.2087 


-0.1611 


0.2861 


0.2187 


-0.7012 


-0.9879 


B(3,4) 


-0.0073 


-0.0742 


-0.0046 


0.0049 


0.0739 


0.0441 


-0.0275 


-0.0068 



Table 9: Predicted coefficients for the OPE and threshold expansion in Eqs. (12.5112.61 
for rii = 3, 4 using the results of Ref [22] . 
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V 


s 




n; = 3 


n t = 4 


ni = 3 


n t = 4 


nj = 3 


n ; = 4 


ni = 3 


m = 4 


n*W(l,0) 


37.5996 


12.0638 


23.3869 


21.2855 


44.7521 


50.5119 


38.2910 


34.1148 


nx(+)(i,i) 


35.4012 


25.9895 








44.8394 


42.4421 








f2*(+)(l,2) 


-11.5083 


-9.8665 








-11.5083 


-9.8665 








nx(.-)(6, i) 


0.6588 


0.4693 


-4.1789 


-3.3437 


-1.4412 


-1.2240 


-9.3411 


-7.7181 


Qx(-)(6, 2) 


-0.6921 


-0.5036 


-0.1225 


-0.2464 


-0.5123 


-0.4476 


-3.6175 


-3.3982 


(6, 3) 


-0.3031 


-0.3166 


0.8217 


0.7098 


0.3593 


0.2710 


1.3931 


0.9956 


Q x <-'(6,4) 


-0.2341 


-0.2016 


0.5560 


0.5083 


0.1412 


0.1334 


3.3030 


3.0905 



Table 10: Predicted coefficients for the OPE and threshold expansion in the oo- variable 
([ZBD for m = 3, 4, using the results of Ref [22]- 
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ni = 3 


ni = 4 


a 2,0 


-19.8034 
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-21.0146 








"2,1 


-12.8916 


-11.0524 








-12.8916 


-11.0524 
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-282.1379 


116.5165 


363.6681 


324.5732 


-1150.5165 


-1030.2552 


2128.4420 


1646.6387 




238.9509 


-120.5650 


-307.8491 


-228.3280 


846.5914 


632.6807 


-2060.9357 


-1931.8120 




82.7196 


-50.8681 


-112.0433 


-100.2002 


372.6561 


334.8105 


-660.5808 


-499.3996 


^3 


-32.4972 


18.2039 


43.6160 


34.0066 


-123.2945 


-97.8701 


260.2617 


231.4261 



Table 11: Induced values for the coefficients in the asymptotic expansion (12. lip for rii = 
3, 4 corresponding to the predicted terms in both the OPE and threshold expansion 
listed in Tables [9] and [KB 





V 






V 




V 




Jl; = 3 
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n v 


ni = 3 


ni = 4 




n t = 3 


n, = 4 


,4(1/2,0) 


-49.6008 


-37.7892 


nc+ 


(1.0) 


111.8000 


22.0766 


02,0 


-4.2693 


-4.4806 


-291.5004 


-106.0299 


1.5961 


-2.0938 


-19.8693 


-18.8806 


A(l/2,1) 


1.7691 


0.8707 


n(+ 


(i,i) 


35.4012 


33.0547 


C(2,l 


-12.8916 


-11.0524 


-13.7341 


-13.4841 


19.8981 


18.6999 


-12.8916 


-11.0524 


.4(1/2,2) 


3.2229 
3.2229 


2.7631 
2.7631 




(1,2) 


-11.5083 
-11.5083 


-9.8665 
-9.8665 


«2,2 








B(3,l) 


-0.9944 


-1.1203 


n<- 


(6, 1) 


1.3333 


1.1538 


/?7,0 


547.6502 


497.1674 


-1.2222 


-1.3309 


0.7307 


0.6148 


-237.6193 


-202.8925 


■8(3,2) 


0.8006 


0.6936 


n<- 


(6,2) 


-0.7378 


-0.6553 


07,1 


585.1540 


541.7054 


0.2451 


0.1830 


-1.2126 


-1.0802 


224.1701 


220.6132 


£(3,3) 


1.2853 


1.1730 


n<- 


(6,3) 


-0.1656 


0.2324 


07,2 


60.6614 


62.7959 


0.3659 


0.3462 


-0.2957 


-0.1236 


-187.7386 


-168.4818 


■8(3,4) 


0.0119 


0.0102 


n(- 


(6,4) 


-0.2357 


-0.2176 


07,3 


-38.4732 


-42.5641 


0.0010 


0.0008 


-0.2330 


-0.2199 


-46.3932 


-45.0972 



Table 12: Predicted coefficients corresponding to the first unknown term in the OPE 
and threshold expansion in the z-variable (second and third columns) and in the in- 
variable (fifth and sixth columns) from the results of Ref [21]. The corresponding coef- 
ficients in the asymptotic expansion (12.111) are also deduced (eight and ninth columns). 
See text. 
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5 Summary and Conclusions 



In this work we have applied a method, previously introduced in Ref. [31], to recon- 
struct all the vacuum polarization functions (12. 2p in the whole complex plane. The 
method uses the first coefficients from the Taylor expansion (12.41) and the nonanalytic 
coefficients of the OPE and threshold expansions (12. 5112. 6ft to resum these three ex- 
pansions simultaneously. The reconstruction is carried out by means of a systematic 
calculation of the asymptotic coefficients Q(n) (12.111) for large n, which control the 
expansion of the vacuum polarization functions in the conformal variable u (12. 3p . This 
reconstruction consists of an approximation to the original function in terms of a finite 
number of polylogarithms (and their derivatives) plus a well-defined polynomial, as 
shown in Eq. (I3.3p . The approximation is controlled by an error function (I3.2p which 
becomes smaller as more information is known. In this regard, we emphasize that 
our present errors are dominated by the unknown terms in the threshold expansion. 
Therefore, a calculation of the coefficients A(l/2, k ^ 0) in the expansion (12. 5 p would 
significantly improve on the precision achieved in the present determination. Alterna- 
tively, knowledge of more terms in the Taylor expansion (I2.4p (i.e. a larger iV*) would 
also help, since currently iV* = 3,4, only. Given that at 0(a 2 s ) nothing less than 30 
terms are known in the vector channel [SJ |9], we may optimistically look forward to an 
improvement on this in the future. 
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Figure 4: Comparison between the prediction for the next Taylor coefficients in the 
pseudoscalar channel and the exact result, as a function of the number of Taylor coef- 
ficients known, N*. See text. 

APPENDIX 



The pseudoscalar channel allows us to check the convergence properties of our 
approximation in the following way. As we did at the end of section 3, we can compute 
the coefficients C p (n) using Eq. (I2.10p by separating the first N* terms in this sum, 
which are exactly known, and using the asymptotic expression (12. lip from the term 
iV* + 1 onwards. By pretending that only one term is known, i.e. N* = 1, one can 
make a prediction for the other coefficients C p (n) with n = 2,3, 4 and then compare 
to their exact value. The result of this exercise is shown in Fig. H] where the black dot 
signifies the exact value known for C p (l) and the arrow leading to the blue dots are 
the predictions for the coefficients C p (2) — C p (4) obtained in this way. Alternatively, 
one may now pretend that two coefficients are known, i.e. iV* = 2, and then predict 
the other two coefficients: this result is represented by the red squares. Finally, if three 
coefficients are assumed to be known, i.e. N* = 3, the prediction for the fourth term is 
the brown triangle. As one can see, there is a very nice convergence as N* is increased, 
already for very low values of iV*. 

The previous exercise was done with the asymptotic expression (I2.1ip including all 
the coefficients a p and (3 P obtained in Table [3] for the pseudoscalar channel. Taking 
C p (l) — C p (3) as exactly known, we may also test the convergence properties of the 
asymptotic expansion f)2.1ip by comparing the result for the coefficient C p (4) and 
the exact result if we now pretend that the last coefficients a^ k ,(k = 0,1,2) and 
/3 p k ,(k = 0,1,2,3) are not knownQ The result is shown in Fig. as a blue square. 
This should be compared to the red triangle, which is the prediction obtained when 

4 The main effect is due to the threshold coefficients a p . The OPE coefficients (3 P play a minor role. 
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Figure 5: Comparison between the prediction for C p (4) and the exact result, as a 
function of the number of terms kept in the asymptotic expansion (12.111) . See text. 

all the coefficients a p and (3 P in Table |3] are included. Again, this shows that the 
asymptotic expression (12. lip converges nicely to the right result. 
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